Blade Displacement Predictions for the Full-Scale UH-60A Airloads Rotor 


Robert T. Biedron 

Senior Research Scientist 

Elizabeth M. Lee-Rausch 

Research Engineer 

NASA Langley Research Center, Hampton, VA 23681 

ABSTRACT 

An unsteady Reynolds-Averaged Navier-Stokes solver for unstructured grids is loosely coupled to a rotorcraft com- 
prehensive code and used to simulate two different test conditions from a wind-tunnel test of a full-scale UH-60A 
rotor. Performance data and sectional airloads from the simulation are compared with corresponding tunnel data to 
assess the level of fidelity of the aerodynamic aspects of the simulation. The focus then turns to a comparison of the 
blade displacements, both rigid (blade root) and elastic. Comparisons of computed root motions are made with data 
from three independent measurement systems. Finally, comparisons are made between computed elastic bending and 
elastic twist, and the corresponding measurements obtained from a photogrammetry system. Overall the correlation 
between computed and measured displacements was good, especially for the root pitch and lag motions and the elastic 
bending deformation. The correlation of root lead-lag motion and elastic twist deformation was less favorable. 


NOTATION 

() tip denotes blade tip 

()oo denotes freestream condition 

a speed of sound [ft/s] 

BD Blade Displacement measurement system 

c section blade chord [ft] 

Cmx hub rolling moment coefficient ( — M 2 g 3 ) 

PooVtip Tin 

CMy hub pitching moment coefficient ( — My 2 3 ) 

^ PooVtip Tin 

Cp rotor power coefficient ( — - p 2 p2 ) 

poo Vfip Tin 

Cp rotor thrust coefficient ( — — 0 ) 

V PooVtip 2 7lR 2J 

Cx rotor propulsive force coefficient ( p y X i nR 2 ) 

CFD Computational Fluid Dynamics 

CSD Computational Structural Dynamics 

f n sectional normal force [lbf/ft] 

f c sectional chord force [lbf/ft] 

m c/4 sectional pitching moment [ft-lbf/ft] 

M 2 C c sectional chord force coefficient ( -r — ) 

6 y \pooCloo 2 c' 

M 2 C m c/4 sectional pitching moment coefficient 


M 2 C n sectional normal force coefficient ( T 9 ) 

n y \pooaoo 2 c J 

M x rotor hub rolling moment [ft-lbf] 

M y rotor hub pitching moment [ft-lbf] 

P rotor power [ft-lbf/s] 
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r radial position [ft] 

R rotor radius [ft] 

T rotor thrust [lbf] 

V velocity [ft/s] 

x chordwise distance from leading edge [ft] 

X rotor propulsive force [lbf] 

a c corrected shaft angle [°] 

a s geometric shaft angle [°] 

u advance ratio ( ^ ) 

v v tip ' 

Poo freestream density [slugs/ft 3 ] 

o solidity 

6 C collective pitch (mean root pitch) [°] 

6\ s sine component of root pitch variation with az- 
imuth [°] 

6\ c cosine component of root pitch variation with az- 

imuth [°] 

1 jf azimuthal position [°] 


INTRODUCTION 

The flight test of the UH-60A airloads rotor in 1993 provided 
an extensive database for the validation of predictive analy- 
sis tools. (Refs. 1,2) Numerous studies have been published 
that compare loosely-coupled Computational Fluid Dynamics 
(CFD) and Computational Structural Dynamics (CSD) pre- 
dictions with the flight test data. The use of CFD has typi- 
cally shown improved comparison of airloads between flight 
test and computation over the lower- fidelity models avail- 
able in CSD comprehensive codes. Despite a high level 
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of fidelity in aerodynamic modeling (including Reynolds- 
Averaged Navier- Stokes (RANS) and Detached-Eddy Simu- 
lations (DES)) and a high level of spatial/temporal resolution, 
the correlation of the predicted and measured sectional air- 
loads has appeared to reach a plateau. (Refs. - 5) Differences 
in sectional airloads between the various simulations are gen- 
erally smaller than the difference between the simulations and 
measured data. A full-scale wind-tunnel test of the same UH- 
60 A rotor used in the 1993 tests was completed in 2010 at 
the National Full-Scale Aerodynamic Complex (NFAC) 40- 
by 80-foot Wind Tunnel. (Ref. 6) The wind-tunnel test pro- 
vided an opportunity to test the airloads rotor in a more con- 
trolled environment without the influence of a tail rotor or 
propulsion system. Also, the wind-tunnel test produced ad- 
ditional data not available from the flight test including new 
measurements and additional flow conditions. The new mea- 
surements include rotor balance forces and moments, oscil- 
latory hub loads, blade displacements and deformations and 
rotor wake measurements. An increased confidence in both 
simulation and measurement can be gained through detailed 
comparisons of this hierarchy of data. Since the completion of 
the NFAC tests numerical simulations at representative condi- 
tions from the test database have begun. Initial simulations 
focused on comparisons with sectional airloads and rotor per- 
formance predictions. (Refs. 7,8) Computations of vortex evo- 
lution have been made using high-resolution structured grids, 
and compared with vortex positions derived from PIV mea- 
surements. (Ref. 9) More recently, computational results fo- 
cusing on blade displacements have been made using OVER- 
FEOW coupled to CAMRAD-II. (Refs. 10, 1 1) 

An unsteady RANS solver for unstructured grids, FUN3D, 
has been used to compute the rotor performance and airloads 
of the UH-60A Airloads Rotor in the NFAC. (Ref. 12) The 
flow solver was loosely coupled to a rotorcraft comprehen- 
sive code, CAMRAD-II, to account for trim and aeroelastic 
deflections. Computations were made for the 1-g level flight 
speed-sweep test conditions with the airloads rotor installed 
on the NFAC Farge Rotor Test Apparatus (FRTA) and in the 
40- by 80-ft wind tunnel to determine the influence of the test 
stand and wind-tunnel walls on the rotor performance and air- 
loads. Detailed comparisons were made between the results 
of the CFD/CSD simulations and the wind tunnel measure- 
ments including rotor loads, sectional loads and blade pres- 
sures. The computed trends in solidity- weighted propulsive 
force and power coefficient matched the experimental trends 
over the range of advance ratios and were comparable to pre- 
viously published results. (Refs. 7, 8, 10, 13-15) Rotor perfor- 
mance and sectional airloads showed little sensitivity to the 
modeling of the wind-tunnel walls, which indicated that the 
rotor shaft-angle correction adequately compensated for the 
wall influence up to an advance ratio of 0.37. Sensitivity of 
the rotor performance and sectional airloads to the modeling 
of the rotor with the FRTA body /hub increased with advance 
ratio. The inclusion of the FRTA in the simulation slightly 
improved the comparison of rotor propulsive force between 
the computation and wind tunnel data at an advance ratio of 
0.30, but did not resolve the difference in the rotor power pre- 


dictions at an advance ratio of 0.37. Despite a more precise 
knowledge of the rotor trim loads and flight condition, the 
level of comparison between the computed and measured sec- 
tional airloads/pressures did not significantly improve over the 
previously published results for the flight test. (Refs. 3, 16) 

The focus of the current paper is to provide detailed com- 
parisons between blade displacements and deflections ob- 
tained from FUN3D/CAMRAD-II simulations and wind tun- 
nel blade displacement measurements from the NFAC test. 
While the emphasis is on blade displacements and deflec- 
tions, comparisons with overall rotor performance measure- 
ments (thrust, power and propulsive force) and sectional air- 
load measurements (normal force, chord force and pitching 
moment) are also presented to provide a more complete as- 
sessment of the fidelity of the numerical simulations for the 
conditions considered in the current paper. 

UH-60A AIRLOADS WIND TUNNEL TEST 

In 2010, a full-scale wind-tunnel test of the UH-60A airloads 
rotor was conducted at the NFAC 40- by 80-foot Wind Tunnel. 
Reference 6 gives a detailed description of the test hardware, 
instrumentation, and data reduction. Figure 1 shows the air- 
loads rotor mounted on the FRTA in the NFAC. The primary 
objective of the test was to acquire a comprehensive set of 
validation-quality measurements on a full-scale rotor at chal- 
lenging conditions outside the conventional flight envelope. A 
secondary objective of the test was to provide data to evaluate 
the similarities and differences between small-scale wind tun- 
nel, full-scale wind tunnel, and full-scale flight test. To meet 
these objectives, the test was separated into six phases: (1) 1-g 
level flight speed sweeps, (2) parametric sweeps, (3) airloads 
flight test simulations, (4) small-scale wind tunnel test simu- 
lations, (5) high advance ratio (slowed rotor) testing, and (6) 
particle image velocimetry (PIV) testing. 



Fig. 1. UH-60A Airloads Rotor mounted on the LRTA in 
the NFAC 40- by 80-foot Wind Tunnel. 

One blade of the airloads rotor was instrumented with 241 
unsteady pressure transducers. Transducers were generally 
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grouped chordwise along nine radial stations. Figure 2 shows 
the radial stations at which test data were collected. A maxi- 
mum of 15 transducers were installed on the upper-surface at 
each radial station, along with a maximum of 15 transducers 
on the lower-surface at each station. As in the flight-test pro- 
gram, the blade pressures were integrated to provide normal 
force, pitching moment, and chord force for each radial sta- 
tion. Although the transducers were clustered near the leading 
edge to better resolve the pressure distribution there, the num- 
ber of points used to integrate the measured data were sparse 
in comparison with the number of points typically used for 
integration of CFD data. Furthermore, not all transducers re- 
mained operational for all of the run conditions of the test 
program. Aerodynamic shear stresses were not measured, so 
the experimental integrated sectional airloads reflect only the 
pressure contribution. 


r/R=0.225 r/R=0.550 r/R=0.775 r/R=0.920 



r/R=0.965 


Fig. 2. UH-60A Airloads Rotor pressure transducer radial 
locations. 

To provide an expanded database for predictive tools, an 
extensive set of measurements were made during the test. 
In addition to the blade instrumented with unsteady pres- 
sure transducers, a second blade was instrumented with strain 
gages. The wind-tunnel test also produced a tremendous 
amount of additional data not available from the flight test in- 
cluding rotor hub forces and moments, oscillatory hub loads, 
blade displacements and deformations, and rotor wake mea- 
surements. The new measurements allow for a more precise 
knowledge of rotor flight and trim conditions. 

The blade displacement (BD) data was obtained from 
multi-camera photogrammetry. (Refs. 10, 17) The photogram- 
metry measurements were based on 48 targets distributed on 
the lower surface of each blade, from radial station r/R = 0.2 
to station r/R = 0.97 at approximately 0.05r/R increments. 
Additional fixed targets in the tunnel ceiling were used for 
calibration purposes. Figure 3 shows the photogrammetry tar- 
gets on the model as installed in the NFAC. 

During data acquisition, three-dimensional target locations 
are projected onto a two-dimensional image plane. Three- 
dimensional target centroid coordinates are then reconstructed 
from multiple two-dimensional images taken from different 
camera positions. Coordinate transforms map the target cen- 
troids from test section coordinates to model coordinates. 
Comparisons with target centroid locations measured on an 
unloaded, undeformed blade then results in estimates of rigid 
and elastic deformations in pitch, flap, and lead-lag (herein 
referred to simply as “lag”). It should be noted that estimates 
of the rigid pitch, flap and lag displacements rely on a data fit 
of the inboard 25 % target locations, along with the assump- 
tion that the blade is effectively rigid at these inboard stations. 



Fig. 3. Lower surface of blade with photogrammetry tar- 
gets used by the Blade Displacement system. 

Elastic deformations that may actually be present at these in- 
board stations may lead to errors in the estimates of the rigid 
pitch, flap and lag angles, which in turn can propagate into the 
estimates of the elastic deformations at other target locations. 
Variations in the measurements from rev-to-rev were noted 
and are included in the data set as standard deviations from 
average values over many revolutions. In general these rev- 
to-rev variations are small compared to blade-to-blade differ- 
ences. As the data reduction for the blade deflections is quite 
complex, the process is still under refinement and the blade 
deflection data used for comparison with the computational 
results should be considered preliminary. 

Blade displacement data from the NFAC tests are cate- 
gorized as primary or secondary. Primary displacement data 
tracks each blade through all four quadrants of rotation. Sec- 
ondary data tracks a different blade in each quadrant. Be- 
cause each blade is different, piecing secondary data together 
to form a complete revolution will result in discontinuities at 
quadrant boundaries. Thus, primary data sets are preferred 
for comparisons with computation. Fortunately, primary data 
sets are available from each of the six test phases mentioned 
earlier. 

Two alternate systems, a crab-arm system and a laser sys- 
tem, also provide independent measurements of the blade-root 
motion. The crab-arm system, which was also employed dur- 
ing the earlier flight-test program, uses a set of three rotary 
variable differential transformers (RVDTs) for each blade. 
The relative motions of the RVDTs are used to obtain blade- 
root motions by solution of three kinematic equations. The 
laser system was comprised of a set of three laser distance 
transducers, one set of three for each blade. Both crab-arm 
and laser data have recently been corrected to account for 


3 




transducer and amplifier drift as well as rotational (centrifu- 
gal force) effects; it is this corrected data that is compared to 
the computational results. 

CFD/CSD METHODOLOGY 

The unstructured-grid Computational Fluid Dynamics (CFD) 
solver used for this study is FUN3D. (Ref. 18) The code 
solves the unsteady RANS equations, with several models 
available for turbulence closure. The solver has a variety of 
mesh-motion options, including rigid, deforming, and over- 
set meshes, and a robust implicit time-advancement scheme. 
For overset meshes, the DiRTlib (Ref. 19) and SUGGAR++ 
(Ref. 20) codes are used to facilitate communication between 
disparate zones in the mesh. For rotorcraft simulations, struc- 
tural dynamics and trim are modeled by coupling with a ro- 
torcraft comprehensive code, a specialized form of Computa- 
tional Structural Dynamics (CSD) code. 

In this paper, an isolated rotor in free air is modeled. An 
unstructured mesh with prismatic elements in the boundary 
layers near blade surfaces, and tetrahedra elsewhere is used to 
discretize the fluid volume surrounding the rotor. The blade 
surface geometry used for grid generation is based on an up- 
dated definition provided by the Sikorsky Aircraft Corpora- 
tion via the UH-60 Airloads Workshop in May 2009. The 
updated geometry includes more spanwise resolution at the 
tip cap and a correction to the spanwise location of the trim 
tab. The near-body blade grid extends away from the blade 
to a cylindrical outer boundary of radius 1.5 c t i p . In the wall 
normal direction, the grid spacing is set such that an average 
normalized coordinate y + is less than one for the first grid cell 
at the wall for the majority of the blade. The maximum spac- 
ing at the blade grid outer boundary is approximately 0. 10 c t i P . 
Each blade grid contains approximately 2.65 million nodes. 
Four blade grids are overset onto a tetrahedral background 
grid to form the computational grid system. 

The background grid is defined in a square box whose 
sides extend 5 R out from the rotor hub. The finest spacing 
in the off-body grid is approximately 0.10 c t i p . This mini- 
mum spacing is maintained, within the constraints of the un- 
structured meshing software, in a cylindrical volume which 
extends I.IR in the blade plane and 0.20R above and 0.20R 
below the blade plane. This cylindrical volume of refinement 
was achieved by using volume sources. The background com- 
ponent grid has approximately 7.0 million nodes. The base- 
line composite grid has a total of approximately 17.6 million 
nodes. A slice through the baseline isolated UH-60 composite 
grid in Fig. 4 shows the spacing characteristics of the blade 
and off-body grids through the computational domain. The 
grids are shown in a vertical plane passing through the rotor 
hub location. The same grid system was used for simulations 
of an isolated UH-60A rotor system in References 3, 12, 16. 

For the results presented here, the one-equation Spalart- 
Allmaras turbulence model is used. (Ref. 21) The turbulence 
model is applied to both the near-body and off-body regions. 
Inviscid fluxes are evaluated with a second-order (in space) 



Fig. 4. Slice through the CFD grid used to model the rotor. 

implementation of Roe’s scheme. (Ref. 22) Time advance- 
ment is carried out with a modified second-order backward- 
difference dual-time scheme which results in reduced leading- 
order truncation error compared to the classical second-order 
backward-difference scheme. (Refs. 23, 24) A time step cor- 
responding to a 1° change in blade azimuth per step is used. 
Time advancement is performed under the guidance of a ‘tem- 
poral error controller’ (Refs. 23,25) in which dual-time subit- 
erations are carried out within each time step until either the 
L 2 norm of the nonlinear residual is reduced to a specified 
fraction of the estimated temporal error, or, a specified max- 
imum number of subiterations have been performed. For the 
current results, the error fraction is 0.05, and the maximum 
number of subiterations allowed is 40. In Reference 3 a de- 
tailed study of the sensitivity of FUN3D results to parameters 
such as grid density, turbulence model, and time step was car- 
ried out for several of the UH-60A flight-test conditions was 
carried out, and the results of that study were used as guide- 
lines for the current simulations. 

The rotorcraft comprehensive code used for this study is 
CAMRAD-II. (Ref. 26) The aerodynamics modules within 
CAMRAD-II are based on lifting-line models using airfoil ta- 
bles, coupled with wake models. Although such aerodynamic 
models can provide reasonable results for many flight con- 
ditions, in some cases the predictions of the airloads can be 
inaccurate because of limitations of the relatively low-order 
aerodynamic modeling. To help overcome these limitations, 
loose coupling to a CFD code may be employed. In the 
loose coupling approach, airloads data from the CFD solver 
and blade motion data from the CSD solver are exchanged 
at relatively infrequent periodic intervals, for example once 
per revolution or more generally once per integer multiple 
of the blade passage frequency. In a typical coupled simu- 
lation, the initial execution of the FUN3D CFD code is car- 
ried out for two complete rotor revolutions, using blade defor- 
mations from a trimmed CAMRAD II solution with unmodi- 
fied lifting-line aerodynamics. In subsequent coupling cycles, 
FUN3D is run for two blade passages between coupling cy- 
cles (i.e., one-half revolution for a four-bladed rotor such as 
the UH-60A). Although coupling could be performed as fre- 
quently as each blade passage, running the solver for the extra 
blade passage helps damp out transients between coupling cy- 
cles. When converged, the loose coupling approach replaces 
the low-order lifting-line aerodynamics of the CSD code with 
the higher-fidelity Navier-Stokes aerodynamics of the CFD 
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code. Loose coupling is appropriate for unsteady problems 
of a periodic nature, such as a rotor fixed within a wind tunnel 
or in steady level flight. A full description of the CFD method 
and coupling strategy can be found in Reference 3. 

Within CAMRAD-II, each blade is modeled as a set of 
nonlinear beam elements. In addition to the structural dy- 
namics modeling, CAMRAD-II offers a sophisticated trim 
capability. For the UH-60A wind-tunnel test conditions pre- 
sented here, a three-degree-of-freedom trim is used, with the 
(solidity-weighted) thrust coefficient, pitching moment and 
rolling moment specified as trim targets within CAMRAD- 
II. The same, widely used CAMRAD-II structural dynamics 
model established in Reference 27 was employed in this study, 
with a couple of modifications that were also used in Refer- 
ences 11, 12. The first change was to increase the number of 
radial aerodynamic panels from 20 to 100. The second change 
was to increase the azimuthal resolution within CAMRAD-II 
from 24 to 360 steps per revolution. As in the standard model, 
seven one-dimensional structural beam elements were used to 
model the rotor blade. 

RESULTS 

During the NFAC tests, data were acquired during all six 
phases of the test campaign to obtain blade displacements and 
deflections with the BD system. Of these, two cases classified 
as primary BD points were numerically simulated and are pre- 
sented in this paper. Run 43 Points 60-63 (herein referred to as 
R43P60) represents a moderate thrust, moderate advance ra- 
tio condition. This run was part of a parametric sweep during 
which the fixed- system hub moments were nominally zero. 
Run 53 Points 20-25 (R53P20) represents a low thrust, high 
advance ratio condition, and was part of the 1-g level flight 
test phase designed to cover representative flight conditions 
of the UH-60A, including representative hub moments. As 
in any test situation, conditions varied slightly from test point 
to test point. Average test conditions for these two data sets 
are summarized in Table 1 . These test points have also been 
examined using OVERLOW/C AMRAD-II: R42P60 in Refer- 
ence 10 and R53P20 in Reference 11. 


Run/Pt 

etc (deg) 


M tip 

Ct / o 

M x ,My (ft-lb) 

R42P60 

0.78 

0.300 

0.649 

0.1002 

-50 My 
44 M x 

R53P20 

-7.99 

0.371 

0.650 

0.0806 

-2,596 My 
-2,840 M x 


Table 1. Average test conditions and trim targets; sign con- 
ventions: a c + aft shaft tilt, My + nose up, M x + right side 
down. 

The CFD/CSD simulations used these average conditions, 
and assume that each blade is identical in terms of mass dis- 
tribution, structural properties, and trim- tab settings. Each ro- 
tor blade in the test is actually different, in ways both known 
(trim-tab setting) and unknown. As will be seen, measure- 
ments show that each blade exhibits different root motions and 


deflections. Note that some of the experimental data presented 
in this paper have been updated since the initial publication in 
References 6, 10. Generally speaking the corrections have 
been small. In the measured data presented here, rev-to-rev 
variations, when available, are shown with error bars. 

Performance and Rotor Trim 

A hierarchy of data is available from the NFAC tests. Perfor- 
mance data reflects integrated values over the entire rotor, and 
therefore comparison between computed and measured per- 
formance data can give an overall assessment of the fidelity of 
the simulation to the test. Details are masked by the integrated 
values however. Table 2 shows the comparison between com- 
putation and measurement for the (solidity-weighted) thrust, 
power, and propulsive force coefficients for the two test con- 
ditions. 




R42P60 

R53P20 

Cy/o 

(NFAC) 

0.1002 

0.0806 

Ct/ g 

(CFD) 

0.1005 

0.0809 

%Diff. 


0.3 

0.4 

Cp/a 

(NFAC) 

0.00382 

0.00845 

Cp/o 

(CFD) 

0.00376 

0.00827 

%Diff. 


-1.6 

-2.1 

C x /o 

(NFAC) 

-0.00218 

0.00956 

C x /o 

(CFD) 

-0.00217 

0.01008 

%Diff. 


-0.5 

5.5 


Table 2. Measured vs. computed performance data. 

The computed power coefficient matches the experimental 
values to within 1.6 and 2.1 percent for R42P60 and R53P20, 
respectively. The propulsive coefficients match to within 0.5 
and 5.5 percent for R42P60 and R53P20, respectively. The 
thrust computed from the CFD solver is within 0.4 percent of 
the measured value for both cases. Note that the target thrust 
(the measured value) is matched essentially exactly within 
the comprehensive code. The small differences in the CFD- 
computed thrust are a measure of losses encountered in a non- 
conservative load transfer process that only transfers three 
(normal force, chord force, and pitching moment) of the six 
force and moment components. 

The trimming process employed here adjusts the collec- 
tive and cyclic pitch components to meet thrust and mo- 
ment targets. The final values from the trim process, ex- 
tracted from CAMRAD-II output, are shown in Table 3. 
The root pitch angle, exclusive of built-in twist and elas- 
tic deformation, is described by 6 = 6 C + Q\ s sin(\\f) + 
Qi c cos(\\f)(+higher hamonics). The collective angles ( G c ) in 
both simulations are slightly larger than the NFAC values, by 
1.12° for R42P60 and by 0.62° for R53P20. 

Sectional Airloads 

As mentioned earlier, a number of pressure transducers be- 
came nonoperational or unreliable during the course of the 
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R42P60 

R53P20 

00 

(NFAC) 

7.88 

13.07 

00 

(CFD) 

9.00 

13.69 

01s 

(NFAC) 

-6.16 

-9.06 

01s 

(CFD) 

-5.83 

-8.47 

01 c 

(NFAC) 

0.41 

1.28 

01 c 

(CFD) 

0.78 

1.00 


Table 3. Measured vs. computed control angles. 

NFAC testing. For Runs 42 and 53, there was not a suffi- 
cient number of operational pressure transducers to provide 
sectional airloads for station r/R = 0.550. Also, for Run 53 
the airloads stations at r/R = 0.675, 0.775, and 0.965 had a 
small number (7-8) of operational transducers on the lower 
surface of the blade. Outboard sectional airloads are gener- 
ally of greater interest owing to higher loading and nonlinear 
phenomenon such as shock wave interactions. Because the 
aforementioned stations from Run 53 do not have sufficient 
transducers to provide viable sectional airloads, a common set 
of outboard stations, r/R = 0.865, 0.920 and 0.990 are exam- 
ined for both test conditions, together with the station closest 
to the blade root, r/R = 0.225. 

The NFAC sectional airloads data is derived solely from 
pressure integration, using a maximum of 30 (and often fewer) 
chordwise transducer locations at each station. In CFD simu- 
lations, sectional airloads are usually integrated over all points 
on the CFD surface at the radial location in question, and the 
integration includes both pressure and skin friction contribu- 
tions. In Reference 16, CFD sectional airloads were obtained 
by integrating only the pressure contribution at locations cor- 
responding to the transducer locations, and only at those loca- 
tions where the transducer was deemed reliable. This pressure 
transducer status information is provided as part of the exper- 
imental data file for each test point. The same second-order, 
mapped-coordinate integration scheme used with the NFAC 
pressures was used with the CFD pressures. As a result, sec- 
tional airloads can be defined consistently in both CFD and 
experiment. Generally speaking this consistent integration 
leads to greatly improved correlation in chord force. The same 
process described in Reference 16 is followed in this paper for 
the computed sectional airloads. 

The measured sectional airloads data shown here reflect a 
correction applied after the initial data release. The prelim- 
inary data used absolute pressures in the integration routine 
rather than relative (P — P static ) pressures that were expected 
by the routine. The resulting corrections were generally small, 
and principally affect the pitching moment. The CFD airloads 
were processed using the relative pressures, and so are con- 
sistent with the corrected NFAC data. The measured sectional 
airloads data are plotted with error bars denoting rev-to-rev 
variations. However, for these test conditions, the error bars 
are almost always smaller than the symbol size used in the 
plots. 

Figure 5 shows the comparison of sectional airloads be- 
tween CFD and NFAC data at radial stations r/R = 0.225, 


r/R = 0.865, r/R = 0.920 and r/R = 0.990 for the test con- 
dition R42P60. Figure 6 shows the corresponding data for 
R53P20. Ordinate ranges for each quantity ( M 2 C n , M 2 C m , 
and M 2 C c ) are held fixed between both test conditions and at 
all three stations, to allow direct comparison of relative mag- 
nitudes. 

The general level of agreement between computation and 
measurement is favorable. At all four stations, the overall 
level of computed normal force, and the variation with az- 
imuth, closely follows the measured normal force. The pitch- 
ing moment is also well captured by the simulations, with the 
exception of r/R = 0.865, where the computation misses the 
magnitude of the nose-down pitching moment on the advanc- 
ing side. Computed chord force, which often exhibits very 
large differences between computation and measurement, fol- 
lows the trend of the measured data quite well, though there 
is a persistent offset from the measured values. As shown in 
Reference 16, this level of agreement in chord force is a direct 
result of the consistent integration of pressures between CFD 
and experiment. 

The high advance ratio case, R53P20, exhibits a pro- 
nounced shock interaction on the advancing side in the out- 
board region of the blade. As the blade rotates through the 
advancing side, the upper surface, normally the suction sur- 
face, and the lower surface change roles, leading to the neg- 
ative loading (negative M 2 C n ) seen in (approximate) azimuth 
range of 90° to 150° at each of the three outboard stations. As 
the shock wave moves from the upper to the lower surface, a 
sharp change in pitching moment is evident near y/ = 75°. 

Overall the level of agreement between simulation and 
measurement suggests that the aerodynamics are generally 
well captured at the outboard blade stations. This is reassur- 
ing, since deflections of the outer part of the blade will depend 
in large part on the aerodynamic loads. Before considering the 
elastic blade deflections, we will examine the blade root mo- 
tions, as these form the basis for defining “rigid blade” mo- 
tions, from which elastic deformations are measured. 

Blade Root Motion 

In this section, the computed blade root motions are extracted 
from sectional geometric data output by FUN3D at each time 
step during a rotorcraft simulation. This geometric data in- 
cludes absolute pitch angle and quarter-chord location in both 
normal and chordwise directions in the blade reference frame. 
The blade reference frame is the coordinate frame used to de- 
fine the unloaded, undeformed blade geometry. This data, in 
conjunction with other geometric data such as built-in blade 
twist, undeformed blade quarter-chord locations, and hinge 
offsets, allows the computed blade motions and deflections to 
be determined. The same type of motion extraction is possible 
not only at the root but at any other radial location at which the 
sectional geometric data is output. FUN3D uses this sectional 
geometric data internally to resolve the sectional airloads into 
normal and chordwise components (as shown in Figure 6 for 
example). Thus the computed blade motion data presented 
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here is obtained in a consistent way with the sectional air- 
loads. 

Comparisons of computed root pitch, root bending and 
root lag motions for R42P60 are shown in Figure 7. Compar- 
ison is made against all three measurement systems. Clearly, 
the measured root motion is different for each blade. In gen- 
eral the measured pitch angles for each blade are within ap- 
proximately 1° of each other, the measured flap angles are 
within approximately 2.5° of each other, and the lag angles 
within approximately 1° to 2.5° of each other. 

The root pitch angle shown in Figure 7 is exclusive of the 
built-in blade twist. The computed root pitch agrees very well 
with the laser data, however the laser data is slightly different 
from the crab-arm and BD pitch data. As observed when dis- 
cussing the trimmed control angles, the computed mean pitch 
angle is roughly 1° higher than that measured by the crab-arm 
and BD systems. 

The computed root flap angle generally falls within the 
bounds of blade-to-blade variation as measured by the crab- 
arm and laser systems. For the root flap angle, the crab-arm 
and laser data are closer to each other than to the BD data, but 
all three systems show similar flap motion. 

The computed lag angles presented here have been in- 
creased by +7°. As discussed in Reference 28 regarding the 
flight-test program, during data acquisition, when the rotor 
hub is at 0° azimuth, the lag angle is +7°. The NFAC tests 
use the same azimuthal reference for data acquisition as the 
flight-test program. In the CFD model, however, the lag angle 
is 0° at the 0° azimuth reference point, thus necessitating the 
+7° increment to lag angle to allow direct comparison with 
the measured value. 

While the predicted pitch and flap angles at the root com- 
pare reasonably well with the measured data, the predicted lag 
angle shows a more significant offset from the data, even ac- 
counting for the previously-mentioned +7° offset. The wave- 
form of the lag response is reasonably well captured by the 
computation, but a +2° to +3° offset is observed between 
computation and measurement, with the largest difference be- 
ing relative to the BD data. This residual lag offset has also 
been observed in References 10, 11. The variation in mean 
lag offset with rotor torque was examined in Reference 1 1 as 
a means of understanding the nature of the offset. There it was 
concluded that one or more of the physical properties of the 
model blade were likely inaccurate, and that errors in torque 
offset, chordwise center of gravity, or tension center were the 
most likely causes leading to the mean lag offset. 

Figure 8 shows the same comparisons of computed and 
measured root motion for the R53P20 (high advance ratio) 
condition. Once again the pitch angle measured by the laser 
system is in best agreement with the computation, though the 
crab-arm and BD data agree better with each other than with 
the laser data. 

Although the excursions in flap angle around the rotor disk 
are larger for R53P20 compared to R42P60, the mean flap 
angle is smaller (3.56° vs. 4.52°), which is consistent with 


the lower thrust level of R53P20. Once again the computed 
flap angles generally lie within the bounds of the blade-to- 
blade variation of the crab-arm and BD measurement, except 
in the range i jf = 300° to i jf = 360°. For this case the laser 
measurements of flap angle are quite different than the crab- 
arm and BD measurements. 

Compared to R42P60, the difference in mean lag angle be- 
tween computation and measurement is smaller, particularly 
compared to the crab-arm and BD data. In fact the com- 
puted root lag angle follows the variation of the BD data 
quite closely, and with the computed lag nearly matching the 
BD data for blade 4. However, the BD data for lag is sus- 
pect. (Ref. 11) 

Blade Deflections 

For both computation and experiment, the evaluation of elas- 
tic blade deformations at a particular azimuth results from the 
subtraction, in the blade reference frame, of a hypothetically 
rigid blade from the instantaneous blade shape. These opera- 
tions are only performed along the blade quarter chord line, so 
that strictly speaking, the “blade deformations” are “deforma- 
tions of the quarter chord”. The instantaneous orientation of 
the hypothetically rigid blade is determined from the instan- 
taneous root pitch, flap, and lag angles discussed in the pre- 
ceding section. Note that because the determination of elastic 
deformations depends on the root motions, errors in the root 
motions will propagate into errors in the elastic deformation. 
This is true for both computation and measurement. 

Consider a plane orthogonal to the axis of rotation, and 
centered on the rotor hub. Flap is a vertical motion out of 
this plane, and lag a motion in this plane. Wherever the blade 
quarter chord differs from the feathering axis, changes in root 
pitch lead to out-of-plane (flap) motion of the hypothetical 
rigid quarter chord in the tip region. On the swept-tip portion 
of the UH-60A blade, this out-of-plane motion due to pitch 
increases linearly with radius, following the increasing dis- 
placement of the quarter-chord line from the feathering axis. 
This out of plane tip motion is accounted for in the hypothet- 
ical rigid blade. 

The elastic component of the flapping motion is referred 
to as elastic bending, while the elastic component of lag mo- 
tion is referred to as chordwise or in-plane bending. Chord- 
wise bending is much smaller in magnitude than elastic bend- 
ing. Whether or not sufficiently accurate measurements of the 
chordwise bending by using the BD system are possible is still 
subject to some debate. No in-plane bending comparisons are 
shown here. 

Figure 9 shows the comparison of predicted and measured 
elastic bending and twist at r/R = 0.97 (just inboard of the tip) 
for the R42P60 conditions. Overall, the elastic bending shows 
very good agreement between the computed data and BD data. 
A phase difference of roughly 20° is observed, but otherwise 
both the amplitude and waveform of the elastic bending are 
very well captured in the computation. The blade-to-blade 
scatter is on the order of 2.5 inches. If the deformed shape of 
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a blade (but not the orientation) was known exactly, an uncer- 
tainty in the root flap angle of 0.5° would lead to a scatter in 
elastic bending of 2.8 inches at the tip of the 322 inch radius 
rotor. 

While the waveform of the elastic twist at r/R = 0.97 is 
reasonably well captured in the simulation, the mean elastic 
twist is more negative, i.e., shows increased washout, com- 
pared to the BD data. This is not too surprising since the mean 
root pitch angle was observed to be greater in the simulation 
than was measured. Since the measured thrust is matched to a 
close approximation by the trim process, a decrease in thrust 
in the outboard sections due to greater washout would need 
to be compensated by an increase in the thrust due to higher 
inboard pitch. 

Figure 10 shows the comparison of predicted and mea- 
sured elastic bending and twist at r/R = 0.97 for the R53P20 
conditions. The same level of agreement between computa- 
tion and measurement is obtained as for the the R42P60 con- 
ditions. The mean tip washout is greater for the higher ad- 
vance ratio conditions, but the offset of the mean washout be- 
tween computation and measurement is roughly the same as 
for R42P60. 

Finally, radial variations in elastic bending and twist are 
compared. Three azimuthal locations are examined: y/ = 0°, 
y/ = 150°, and y/ = 255°. The last two correspond, to the 
(approximate) azimuthal location where the most negative 
(downward) elastic tip bending occurs, and where the most 
positive (upward) tip bending occurs, respectively. 

Figures 1 1 and 12 show a comparison of radial variation in 
elastic bending and twist for the R42P60 and R53P20 test con- 
ditions, respectively. The computed elastic bending is in rea- 
sonably good agreement with the measured bending for both 
test conditions. As would be expected, the radial variation 
of elastic bending primarily consists of a first beam bending 
mode, although higher bending modes are clearly present. 

The correlation between computation and measurement of 
elastic twist as a function of radius is considerably less favor- 
able than that for elastic bending, particularly at y/ = 150°. 
As seen in Figures 7 and 8, yr = 150° is in the region where 
the difference in root pitch between computation and BD data 
is the largest. The elastic twist data at yr = 0° for blade 1 
under the R42P60 condition certainly seems suspect because 
the values actually increase toward the root. The radial dis- 
tributions show that the differences in tip elastic twist seen in 
Figure 9, are not just localized to the tip, but are present across 
the entire blade. 

As mentioned earlier, errors in the rigid motions at the 
root will propagate into errors in the elastic deflections. Al- 
though not carried out here, estimated corrections can be ap- 
plied to computed deflections by modifying the hypothetical 
rigid blade notion with offsets arising from the differences be- 
tween measured and computed rigid root motions, as done in 
Reference 11. 


CONCLUDING REMARKS 

The wide range of data being made available from the full- 
scale tests of a UH-60A rotor in the NFAC are proving to be 
extremely valuable in evaluating computational models. The 
range of data available, from global measures such a rotor 
power and propulsive force, to finer-grain measures such as 
sectional airloads and blade deflections, down to the level of 
unsteady pressures, provides a useful hierarchy for computa- 
tional validation. 

In this paper we explored several levels of this data hier- 
archy, with emphasis placed on correlation of computed and 
measured blade displacements. Both rigid motion at the root 
and elastic deflections along the blade were examined. Two 
test conditions were modeled, one a moderate thrust, moder- 
ate advance ratio case with nulled pitch and roll moments, and 
a second representative of 1-g level flight at a high advance ra- 
tio. 

In terms of overall performance measures, the measured 
power coefficients for the high advance ratio case were com- 
puted to within approximately two percent, and the computed 
propulsive force to within approximately five percent. For the 
moderate advance ratio case the computed results for propul- 
sive power were considerably better. 

Sectional airloads were extracted from the computation in 
a manner consistent with the tunnel measurements, using only 
the pressure contribution, and only at locations corresponding 
to operational pressure transducers. The correlation between 
computation and measurement was generally very good. 

Blade root motions were compared with data from 
three independent measurement systems, one of which, the 
photogrammetry-based blade displacement system, was used 
to obtain displacements along the blades as well as at the root. 
The cases chosen for this paper were ones for which BD data 
were taken for each of the four blades over the entire rotor 
disk. Considerable blade-to-blade variation in the measured 
data is observed, particularly in flap and lag. Typically the 
data from the crab-arm and BD systems agreed with each 
other better than the laser system. Root pitch and flap were 
reasonably captured by the simulations, though - and this is 
most readily apparent for the higher advance ratio - there is 
a phase shift of approximately 20° between the BD data and 
computed results for both pitch and flap. The cause of this 
shift is not yet understood. 

Elastic bending and twist are determined by subtracting 
from the instantaneous blade shape a hypothetically rigid 
blade that moves with the blade-root motion. At the tip, where 
elastic deformations are largest, the elastic bending predicted 
by the computation is in excellent agreement with the data, 
apart from the 20° shift also seen in the root motion. Elastic 
twist on the other hand, is consistently underpredicted in the 
simulations. 

Because the elastic motions are derived motions relying on 
the accuracy of the root motion, future work perhaps should 
examine the absolute position and orientation of the blade, 
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rather than introduce the artifice of the elastic deflections rel- 
ative to a hypothetically rigid blade. 
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Fig. 7. Comparison of predicted and measured blade root motion, R42P60 (/i = 0.300, Cj/cr = 0.1002, a c = 0.78°). 
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Fig. 11. Comparison of radial variation in elastic bending and elastic twist at selected azimuthal locations, R42P60 
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Fig. 12. Comparison of radial variation in elastic bending and twist at selected azimuthal locations, R53P20 
(ji = 0.371, C T /o = 0.0806, a c = -7.99°). 
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